Superconducting symmetry of three-dimensional t-J model on simple cubic lattice 
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Motivated by the finding of nearly isotropic superconductivity in (Ba, K)Fe2As2, we use renormal- 
ized mean field theory to investigated the t-J model on three-dimensional simple cubic lattice. A 
tunable anisotropic parameter is introduced to dictate the coupling on z direction. The symmetry of 
the superconducting order is studied in detail. Calculation shows that for the isotropic case, pairing 
parameters on the three perpendicular directions have |7r phase shift to each other. However, when 
the interaction on z direction is suppressed, the corresponding amplitude of the pairing parameter 
decreases rapidly, furthermore, two-dimensional d-wave state pairing is favored when the anisotropic 
rate less than 0.75. 

PACS numbers: 74.20.Rp, 74.20.Mn 



I. INTRODUCTION 



The key feature of copper oxides is the layered struc- 
ture and led to speculation that reduced dimensional- 
ity is a necessary prerequisite for superconductivity at 
temperatures above iOK [l|. Although two dimensional 
(2D) models, such as t-J model or Hubbard model, 
have captured essential of superconductivity, and suc- 
cessfully explained properties of the un-doped insula- 
tor and occurrence of gap in superconductors, 2D mod- 
els alone can not describe and explain all observations 
of experiments 0, H, H, H, 0, 0|- Despite how large 
the ratio between out-of-plane and in-plane resistivity 
is, at the phase transition temperature Tc, both resis- 
tivities drop to zero simultaneously, this indicates the 
phase transition is of three dimensional 0], meanwhile, 
the observed antiferromagnetism is definitely a 3D phe- 
nomenon Some works provided evidence that super- 
conductivity in the infinite-layer compounds ACu02 
is of three-dimensional nature which do not contain a 
charge-reservoir block and the distance from one unit 
cell to the next is the shortest among all the cuprates. 
Experimental data 0, H, H, H, H, 0| shows that the de- 
crease of doping concentration is accompanied by a raise 
of anisotropy which is defined as the ratio of the cor- 
relation lengths parallel and perpendicular to the Cu02 
plane in cuprates superconductivity. For materials Y123, 
Y124 0, m and HgBa2Cao.86Sro.i4Cu206-5 [3, co- 
herence length and the anisotropy ratio imply that they 
are anisotropic 3D superconductors, these phenomena 
are supported b y g ood 3D scaHng analysis [13 • Re- 
cently observed [l3l . [3, superconductivity in iron- 
arsenic-based compounds has attracted many atten- 
tion. Our motivation of investigating three dimensional 
isotropic superconductivity come directly from the mea- 
surements of the electrical resistivity in single crystals of 
(Ba, K)Fe2As2 in a magnetic field up to GOT [l^, Yuan et 
al found that the superconducting properties are in fact 
quite isotropic, appear more three dimensional than that 
of the copper oxides. Their results indicates that reduced 
dimensionality in these compounds is not necessarily a 



prerequisite for high temperature superconductivity. 

3D anisotropic t-J model has already been studied be- 
fore By using mean-field Hamiltonian and 
carrying out ex pan sion of free-energy, two main results 
were obtained [17| : One is that transition temperature 
decreases weakly with both increasing of 3D coupling 
strength and doping concentration, the other one is that 
in all cases d-wave pairing ansatz has the lowest energy. 
However, in simple cubic lattice (SCL), each site has six 
nearest neighbors (nn) settled in three perpendicular di- 
rections, no direction is special, if superconducting be- 
havior is possible, its symmetry can not be conventional 
d-wave. 

With the help of renormalized mean- field theory [H 
|2l|(RMFT), We found that for isotropic case, pairing 
parameters on the three perpendicular directions have 
|7r phase shift to each other. While as interaction in z 
direction is suppressed, corresponding amplitude of pair- 
ing parameter drops quickly from infinite value to zero. 
By tuning coupling integral in z direction, our calcula- 
tion shows that superconducting symmetry are functions 
of the anisotropic parameter and doping concentration. 
Moreover, as the anisotropic parameter decreases from 1 
to 0.75, symmetry of pairing parameters change from ^tt 
of 3D to d-wave of 2D. This may give some understand- 
ing of 3D-2D crossover. 



II. FORMULATION 

In SCL, anisotropic t-J model can be written as H = 
PdHtPd + Hs with 

(nn)cr (nn_|_)cr 

= JJ2S^■S,+ J77 ■ 5„ (1) 

(nn) (nnj^) 

where Pd = ]^(1 — rii^riii) is the Gutzwiller projection 



operator [22, |23| which removes totally the doubly occu- 
pied states, t and J are the electron hoping interaction 
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and antiferromagnetic exchange interaction, respectively. 

is to create an electron with spin a at site i, and Si 
is a spin operator. Summation (nn) runs over all nn 
in xy plane, while summation {nn\^ runs over all nn 
in direction z which is perpendicular to xy plane. For 
convenience, all anisotropic parameters are put into z di- 
rection terms, A and 77 are anisotropic parameters with 
range [0,1], A = 77 = 1 corresponding to the isotropic 
case . 

In RMFT the wavcfunction of the Hamiltonian is 
assumed to be the projected state IvE*) — Pd\'^ bcs) , 

strained in the reduced Brillouin zone, and the two co- 
efficients satisfy \uk\'^ + I'-'fcP ~ 1- The projection op- 
erator can be taken into account by a set of rcnormal- 
ized factors [13, HI] defined as (c|g.Cjo-) ~ gt{c\^Cjc)o, 
{Si ■ Sj) « gs{Si ■ Sj)o, where ()o denotes expectation 
value of unprojected state ^'bcs, and () denotes expec- 
tation value of physical state Then one has (H) — 
{H')q — {gtHt + gsHs)o- In homogenous case the renor- 
malized factors [20, [23 take the form of = 2(5/(1 + S) 
and = 4/(1 + 6)'^. Considering even-parity case in 
which u_j:v* ~ = Uj^v~ and = the expecta- 

tion value of the effective hamiltonian has the same form 
as that of 2D 

{H')o = 2gtY,e^,\vj:\^+N-' 

k 

k,k' 

where Ns is the total number of sites and 




FIG. 1: Parameters amplitude as functions of doping concen- 
tration S. ^ = denotes amplitude of particIe-hoIe parame- 
ter, A — \ Ai\ denotes pairing parameter, and As is supercon- 
ducting parameter defined as ptA. 



77^2 cos k^, 
-)o. These 



rjA^coskz, ^ £k ~ ixcoskx - iyCosky - 

h = - m)/(|.9s^), and = + N-\^ 
gap equations should be solved simultaneously with dop- 
ing concentration 5 = ^r^J2k^k/-^k- ^ftcr iterative 
self-consistent solving, for a set of given 6, A. 77 and /2 one 
can obtain all those particle-particle and particlc-holc pa- 
rameters simultaneously. Superconductivity symmetry is 
determined by the phase shift of different pairing param- 
eters At- and the superconductivity parameter A^t f20| 
equals gtAT. 



Vk 

V 

'k 



2 (cos kx -f cos ky + A cos fc^), 
2 (cos kx -+ cos ky -f 77 cos kz). 



(3) 



In order to investigate superconducting property, 
one should introduce two mean-field parameters 
such as particlc-particlc(pairing) parameter A,- = 

(''it'^i+Ti ~ '^Ix'^I+tt)" ^^'^ particle-hole parameters = 
X]cr(cLci+T,(T>o- By minimizing the quantity W ^ {H' - 
MSiiT c|crCicr)o with rcspcct to uj: and f^, where /i is de- 
noted as chemical potential, one gets the coupled gap 
equations 



A, = Ns^Y.^oskrAj:/Ej:, (4) 

k 
k 

where r indicates the three perpendicular directions 



x,y,z, Ef;. = J^l 



A,r COS k,r 



Ay cos ky 



III. SYMMETRY OF SUPERCONDUCTIVITY 
FOR ISOTROPIC AND ANISOTROPIC CASES 

For isotropic SCL, 77 = A = 1. In half-filled case, 
p,, £k = 0, there is a trivial solution with A^ = 
corresponding to projected fermi-liquid. changes its 
sign at the surface and the average energy of per site 
is = -^E{£,jj,<o7fc_fe/ ~ -0.505J. For non- 
trivial solution, by using Eq.®, the energy per site 
can be written as lj = ^fsst/X^fc^fci here relations 
l"fel^ = ^(1 + S,f:/Ej:) and Ujtwl = have been used. 

By assuming E^ = c(cosfca;^ -|- cosfcy^ -|- cosfc^^)^/^ and 
substituting it into gap equations, one can get c = 
3^Sfc(cosfca;^ + COS ky^ + cosfcz^)^/^ ~ 0.398 and the 
energy of per site is w = — |c^J w — 0.712J, which is 
lower than the energy of the projected fermi liquid state 
and is more favored and stable. In the non-trivial case 
parameters should satisfy following equations simultane- 
ously 



(A, 



,A;^+h.c) + 2Cr,Cr, 



(6) 
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It has SU(2) degeneracy, the most important solution is 




= 0.229, 



|Ar| = = 0.324, 

= |A|expi6', 
Ay = |A|expi(6' + 2/37r), 
A;, = |A|expi(6' + 4/37r), (7) 

It clearly shows that the phase shift of different A,- is 
^ . By changing the sign of S^y and taking the phase dif- 
ference of any two pairing parameters as |^7r|, one can 
obtain another solution and if one sets one or two of the 
three ^ as zero other solutions can also be obtained. All 
these solutions have the same energy. Among these ener- 
getically degenerated states the |7r symmetric state has 
the best kinetic energy {Ht)o- Upon doping degeneracy 
will be lifted, superconducting state favors the best ki- 
netic energy state, which is the ^ symmetry state. This 
is also the reason why we call this solution as the most 
important one. 

Hoping integral t is used as energy unit, and t/ J = 3 is 
taking in order to be consistent with the supercxchange 
relation of J = 4t^/?7 in the large Hubbard U limit. For 
isotropic SCL case, self-consistent parameters as func- 
tions of doping concentration are shown in Fig[TJ Am- 
plitude of all At are the same which is denoted as A in 
the figure. Every is real and has the same value of 
^. One can see from Fig[T]that with doping increasing, 
amplitude of pairing parameters decreases, while super- 
conducting parameter A^ varies along a non-monotonic 
curve. These properties are similar to that of 2D square 
lattice. The most interesting result is that each At has 
imaginary part, 9r is used to denote phase of A,-, the 
pairing parameters have phase shift to each other just 
as that of half- filled case. 

When the interaction in z direction is suppressed, am- 
plitude of the corresponding parameters will deviate from 
those of xy plane. In order to make the situation more 
simpler, we set 77 = A. For j] = 0.9, the doping depen- 
dent parameters are presented in Fig. [2l Fig. ^a.) shows 
that 77 affects amplitudes of both paring parameter and 
particle-hole parameter, with increasing S all pairing pa- 
rameters decrease with jA^j < jA^i-^j^j. Anisotropy also 
affects the symmetry of the A^-. Fig. [2jb) demonstrates 
that at half-filled point 9^ y = 9^ — 9y < |7r, and O^^y 
decreases with increasing S. Accompanied by decrease 
of lA^I, 9x,y approaches to tt. For 77 = 0.8, as shown in 
Fig. [UJa) with S increasing jA^I drops more rapidly than 
|Aa;^j,| and vanishes at S — 0.1. Symmetry of pairing pa- 
rameters are shown in FiglSjb), 9x,y decreases from about 
I.OStt to tt at (5 = 0.1. For 5 > 0.1, system apparently 
behaves as 2D with superconducting order being d^^-y^ 
symmetry. By compare above two anisotropic cases one 
can reasonably expect that at a given anisotropic param- 
eter, system will behave as 2D in all doping level. 
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FIG. 2: Left picture shows mean-field parameters amplitude 
as function of S for 77 = A = 0.9. Right picture shows how the 
phases 9t varied as function of 5. 
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FIG. 3: Left picture shows mean-field parameters amplitude 
as function of S for rj = X — 0.8. Right picture shows how the 
phases 9t varied as function of 5. 



This property can be demonstrated clearly in half- 
filled case. For anisotropic half-filling, Ej: = 
•y/ cos^ kx + cos^ ky + c^rf cos^ where = + 
|A|f , c§ = ^? + |A|^, i represents x 01 y direction. For a 
given rj one can obtain the phase difference 9x^y and the 
value of ^ for the best kinetic energy state. As C3 = c, 
it reduces to the isotropic case. As C3 approaches to 0, 
degree of anisotropy is very large and the system turns 
to a quasi-2D one. From Fig. [4] one can see that by 
decreasing 77 from 1, ^ decreases quickly and reaches 
zero at about 77 = 0.75, simutanously the phase differ- 
ence 9x^y increases from |7r to tt. These results indicate 
that as the anisotropic parameter decreases to 0.75, A^ 
vanishes, and system loses its 3D character. 
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FIG. 4: For half-filled anisotropic case, with 77 decreasing, 
amplitude of C3 = [Ajp + Cz drops quickly. Phase difference 
9x^y as function of anisotropic coefficient is presented in right 
picture. 

IV. SUMMARY 

Experiment shows that 122-type ternary iron arsenides 
possess three-dimension properties [l6j , although the mi- 
croscopic nature of superconductivity in iron-based com- 



pound is not clearly at present, and one band model is not 
enough to describe them, we investigated isotropic and 
anisotropic t-J model on simple cubic lattice to show 
the superconductivity symmetry from mean-field point 
of view. For isotropic three-dimensional t-J model, su- 
perconductivity ground state is not conventional d-wave, 
phase shift of each pairing parameter is exactly |7r. For 
anisotropic cases three-dimensional character is not so 
obviously, adding a small anisotropic interaction on z di- 
rection will induces a great anisotropy in its correspond- 
ing mean-field paring parameter and raise serious insta- 
bility of previous 3D superconducting symmetry. We 
found that pairing parameter depends strongly on 
the anisotropic parameter, as anisotropic parameter de- 
crease to 0.75, system appears 2D behavior. From this 
discussion one can see that 3D character superconductor 
is sensitive to amplitude of couplings. 
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